Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.

Evaluating functional network inference using simulations of complex biological systems.

MOTIVATION: Although many network inference algorithms have been presented in the bioinformatics literature, no suitable approach has been formulated for evaluating their effectiveness at recovering models of complex biological systems from limited data. To overcome this limitation, we propose an approach to evaluate network inference algorithms according to their ability to recover a complex functional network from biologically reasonable simulated data. RESULTS: We designed a simulator to generate data representing a complex biological system at multiple levels of organization: behaviour, neural anatomy, brain electrophysiology, and gene expression of songbirds. About 90% of the simulated variables are unregulated by other variables in the system and are included simply as distracters. We sampled the simulated data at intervals as one would sample from a biological system in practice, and then used the sampled data to evaluate the effectiveness of an algorithm we developed for functional network inference. We found that our algorithm is highly effective at recovering the functional network structure of the simulated system-including the irrelevance of unregulated variables-from sampled data alone. To assess the reproducibility of these results, we tested our inference algorithm on 50 separately simulated sets of data and it consistently recovered almost perfectly the complex functional network structure underlying the simulated data. To our knowledge, this is the first approach for evaluating the effectiveness of functional network inference algorithms at recovering models from limited data. Our simulation approach also enables researchers a priori to design experiments and data-collection protocols that are amenable to functional network inference.

Sea Urchin Body Plan Development and Evolution: An Integrative Transcriptomic Approach

My dissertation work integrates comparative transcriptomics and functional analyses to investigate gene expression changes underlying two significant aspects of sea urchin evolution and development: the dramatic developmental changes associated with an ecologically significant shift in life history strategy and the development of the unusual radial body plan of adult sea urchins.

In Chapter 2, I investigate evolutionary changes in gene expression underlying the switch from feeding (planktotrophic) to nonfeeding (lecithotrophic) development in sea urchins. In order to identify these changes, I used Illumina RNA-seq to measure expression dynamics across 7 developmental stages in three sea urchin species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph Heliocidaris tuberculata, and an outgroup planktotroph Lytechinus variegatus. My analyses draw on a well-characterized developmental gene regulatory network (GRN) in sea urchins to understand how the ancestral planktotrophic developmental program was altered during the evolution of lecithotrophic development. My results suggest that changes in gene expression profiles occurred more frequently across the transcriptome during the evolution of lecithotrophy than during the persistence of planktotrophy. These changes were even more pronounced within the GRN than across the transcriptome as a whole, and occurred in each network territory (skeletogenic, endomesoderm and ectoderm). I found evidence for both conservation and divergence of regulatory interactions in the network, as well as significant changes in the expression of genes with known roles in larval skeletogenesis, which is dramatically altered in lecithotrophs. I further explored network dynamics between species using coexpression analyses, which allowed me to identify novel players likely involved in sea urchin neurogenesis and endoderm patterning.

In Chapter 3, I investigate developmental changes in gene expression underlying radial body plan development and metamorphosis in H. erythrogramma. Using Illumina RNA-seq, I measured gene expression profiles across larval, metamorphic, and post-metamorphic life cycle phases. My results present a high-resolution view of gene expression dynamics during the complex transition from pre- to post-metamorphic development and suggest that distinct sets of regulatory and effector proteins are used during different life history phases.

Collectively, my investigations provide an important foundation for future, empirical studies to investigate the functional role of gene expression change in the evolution of developmental differences between species and also for the generation of the unusual radial body plan of sea urchins.

Engineering Highly-functional, Self-regenerative Skeletal Muscle Tissues with Enhanced Vascularization and Survival in Vivo

Tissue engineering of biomimetic skeletal muscle may lead to development of new therapies for myogenic repair and generation of improved in vitro models for studies of muscle function, regeneration, and disease. For the optimal therapeutic and in vitro results, engineered muscle should recreate the force-generating and regenerative capacities of native muscle, enabled respectively by its two main cellular constituents, the mature myofibers and satellite cells (SCs). Still, after 20 years of research, engineered muscle tissues fall short of mimicking contractile function and self-repair capacity of native skeletal muscle. To overcome this limitation, we set the thesis goals to: 1) generate a highly functional, self-regenerative engineered skeletal muscle and 2) explore mechanisms governing its formation and regeneration in vitro and survival and vascularization in vivo.

By studying myogenic progenitors isolated from neonatal rats, we first discovered advantages of using an adherent cell fraction for engineering of skeletal muscles with robust structure and function and the formation of a SC pool. Specifically, when synergized with dynamic culture conditions, the use of adherent cells yielded muscle constructs capable of replicating the contractile output of native neonatal muscle, generating >40 mN/mm2 of specific force. Moreover, tissue structure and cellular heterogeneity of engineered muscle constructs closely resembled those of native muscle, consisting of aligned, striated myofibers embedded in a matrix of basal lamina proteins and SCs that resided in native-like niches. Importantly, we identified rapid formation of myofibers early during engineered muscle culture as a critical condition leading to SC homing and conversion to a quiescent, non-proliferative state. The SCs retained natural regenerative capacity and activated, proliferated, and differentiated to rebuild damaged myofibers and recover contractile function within 10 days after the muscle was injured by cardiotoxin (CTX). The resulting regenerative response was directly dependent on the abundance of SCs in the engineered muscle that we varied by expanding starting cell population under different levels of basic fibroblast growth factor (bFGF), an inhibitor of myogenic differentiation. Using a dorsal skinfold window chamber model in nude mice, we further demonstrated that within 2 weeks after implantation, initially avascular engineered muscle underwent robust vascularization and perfusion and exhibited improved structure and contractile function beyond what was achievable in vitro.

To enhance translational value of our approach, we transitioned to use of adult rat myogenic cells, but found that despite similar function to that of neonatal constructs, adult-derived muscle lacked regenerative capacity. Using a novel platform for live monitoring of calcium transients during construct culture, we rapidly screened for potential enhancers of regeneration to establish that many known pro-regenerative soluble factors were ineffective in stimulating in vitro engineered muscle recovery from CTX injury. This led us to introduce bone marrow-derived macrophages (BMDMs), an established non-myogenic contributor to muscle repair, to the adult-derived constructs and to demonstrate remarkable recovery of force generation (>80%) and muscle mass (>70%) following CTX injury. Mechanistically, while similar patterns of early SC activation and proliferation upon injury were observed in engineered muscles with and without BMDMs, a significant decrease in injury-induced apoptosis occurred only in the presence of BMDMs. The importance of preventing apoptosis was further demonstrated by showing that application of caspase inhibitor (Q-VD-OPh) yielded myofiber regrowth and functional recovery post-injury. Gene expression analysis suggested muscle-secreted tumor necrosis factor-α (TNFα) as a potential inducer of apoptosis as common for muscle degeneration in diseases and aging in vivo. Finally, we showed that BMDM incorporation in engineered muscle enhanced its growth, angiogenesis, and function following implantation in the dorsal window chambers in nude mice.

In summary, this thesis describes novel strategies to engineer highly contractile and regenerative skeletal muscle tissues starting from neonatal or adult rat myogenic cells. We find that age-dependent differences of myogenic cells distinctly affect the self-repair capacity but not contractile function of engineered muscle. Adult, but not neonatal, myogenic progenitors appear to require co-culture with other cells, such as bone marrow-derived macrophages, to allow robust muscle regeneration in vitro and rapid vascularization in vivo. Regarding the established roles of immune system cells in the repair of various muscle and non-muscle tissues, we expect that our work will stimulate the future applications of immune cells as pro-regenerative or anti-inflammatory constituents of engineered tissue grafts. Furthermore, we expect that rodent studies in this thesis will inspire successful engineering of biomimetic human muscle tissues for use in regenerative therapy and drug discovery applications.

Changes in Brain Resting-state Functional Connectivity Associated with Peripheral Nerve Block: A Pilot Study.

BACKGROUND: Limited information exists on the effects of temporary functional deafferentation (TFD) on brain activity after peripheral nerve block (PNB) in healthy humans. Increasingly, resting-state functional connectivity (RSFC) is being used to study brain activity and organization. The purpose of this study was to test the hypothesis that TFD through PNB will influence changes in RSFC plasticity in central sensorimotor functional brain networks in healthy human participants. METHODS: The authors achieved TFD using a supraclavicular PNB model with 10 healthy human participants undergoing functional connectivity magnetic resonance imaging before PNB, during active PNB, and during PNB recovery. RSFC differences among study conditions were determined by multiple-comparison-corrected (false discovery rate-corrected P value less than 0.05) random-effects, between-condition, and seed-to-voxel analyses using the left and right manual motor regions. RESULTS: The results of this pilot study demonstrated disruption of interhemispheric left-to-right manual motor region RSFC (e.g., mean Fisher-transformed z [effect size] at pre-PNB 1.05 vs. 0.55 during PNB) but preservation of intrahemispheric RSFC of these regions during PNB. Additionally, there was increased RSFC between the left motor region of interest (PNB-affected area) and bilateral higher order visual cortex regions after clinical PNB resolution (e.g., Fisher z between left motor region of interest and right and left lingual gyrus regions during PNB, -0.1 and -0.6 vs. 0.22 and 0.18 after PNB resolution, respectively). CONCLUSIONS: This pilot study provides evidence that PNB has features consistent with other models of deafferentation, making it a potentially useful approach to investigate brain plasticity. The findings provide insight into RSFC of sensorimotor functional brain networks during PNB and PNB recovery and support modulation of the sensory-motor integration feedback loop as a mechanism for explaining the behavioral correlates of peripherally induced TFD through PNB.